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This paper studies baryons and baryonic matter in the combined large N c and heavy quark 
mass limits of QCD in 1+1 dimension. This regime allows for a simple physical picture, which is 
computationally tractable. In this non-relativistic limit, baryons are composed of N c quarks that 
interact, at leading order in N c , through a color Coulomb potential. Using variational techniques, 
single baryon masses and interaction energies of low-density baryon crystal are calculated. The single 
baryon calculations are very accurate and can be used as a cross-check to the general numerical 
approach based on Hamiltonian methods for baryons and baryonic matter in 1+1 dimension for 
arbitrary quark mass recently proposed by Bringoltz, which is based on a lattice in a finite box. It is 
noteworthy that the Bringoltz method differs from a previous approach developed by Salcedo, et al. 
in its treatment of an effect due to the finite box size with periodic boundary conditions - namely, the 
existence of gauge configurations which wind around the box. As these are finite volume effects, one 
might expect them to be small for boxes that are large enough so that the baryon density approaches 
zero to high accuracy at the edges of the box. However, the effects of these windings appear to be 
quite large even in such boxes. The large mass infinite volume calculations performed here are 
consistent with the results of numerical calculations using the Bringoltz approach and confirm the 
importance of gauge fields that wrap around for finite systems. The calculation of the baryon crystal 
interaction energy requires the assumption that at low-densities the ground state is composed of 
individual baryons, each in color-singlet states, with wave functions orthogonal to one another. 
While plausible, this assumption is somewhat ad hoc in that one can construct configurations in 
which the entire state is color-singlet but cannot be broken into individual color-singlet baryons. 
However, the interaction energy of low-density baryon crystals calculated with this assumption is 
consistent with numerical results based on Bringoltz's approach. This suggests that the assumptions 
regarding the color structure is justified in 1+1 dimension. This is useful, as it supports a similar 
assumption that was made in 3+1 dimensions, where there is no alternative means of calculation. 

PACS numbers: 11.15.Pg, 12.39.Hg, 14.20.-c, 21.65.Mn 



I. INTRODUCTION 



The problem of understanding cold nuclear matter 
from QCD first principles has remained quite challeng- 
ing due to the intrinsic non-pcrturbative nature of the 
problem [l], Q ■ While many non-perturbative problems 
such as the finite temperature problem can be studied 
using lattice QCD techniques, the problem of finite den- 
sity QCD is marred by the presence of the fermion sign 
problem [l|, ; only low densities can be studied and 
that too only at large temperatures Q as all available 
methods are only valid for small values of ^. As such 
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it is useful to find certain limiting cases for the study of 
cold nuclear matter, in which the problem of finite den- 
sity QCD becomes somewhat more tractable. 

Large N c QCD is a limit where the physics simpli- 
fies considerably 0- Only planar diagrams contribute 
at leading order in N c while quark loops effects are sup- 
pressed. Despite these simplifications, the problem of 
large N c QCD has remained unsolved. However, the 
problem in 1+1 dimension is more tractable [1, here 
the physics simplifies considerably because gluons, due to 
the absence of magnetic fields in one spatial dimension, 
do not propagate. Additionally, with appropriate choice 
of gauge (either axial or Coulomb), gluon-gluon interac- 
tions are absent and hence the theory looks very much 
like QED, in which color charges (i.e. quarks) interact 
via a Coulomb potential. The self-energy corrections to 
the quark propagator come from an infinite number of 
ladder diagrams, which remarkably can be summed up 
exactly in spite of the need to solve a bootstrap equation. 
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Furthermore, the meson spectrum can be solved - albeit 
only numerically via a Bethe-Salpeter equation; a linear 
"Rcggc trajectory" (with no continuum) is observed. 

Unlike the meson problem in 1+1 large iV c QCD, the 
baryon problem is highly non-trivial. Man y a ttempts 
have been made to study the problem: Ref. [l(| studies 
the problem in the chiral limit, where baryons assume the 
form of a Sine-Gordon soliton. Ref. [ll| invokes Eguchi- 
Kawai reduction [l2| to argue controversially jl3l | that 
1+1 large N c QCD is independent of baryon chemical po- 
tential. Ref. [l4| formulates the large N c baryon problem 
for arbitrary quark masses at infinite volume. Numerical 
solutions of the problem require a lattice treatment in 
a finite-sized box. The most exhaustive work was done 
recently by Bringoltz in Ref. 1131 . It is based on the axial 
gauge representation of QCD |15| and the coherent states 



formalism [16|. Ref. [13( solves both the single baryon 



and baryon crystal problems in a finite sized box for ar- 
bitrary quark masses. The principal difference between 
the approaches of refs. [l~3| and [l4[ is the treatment of 
the effects of finite box sizes. In Ref. [l3| the infinite 
volume Hamiltonian is derived and simply used on the 
finite box by imposing boundary conditions. In contrast, 
the idea underlying Ref. [13| is that the quarks making 
up a baryon in a finite box interact via a color-Coulomb 
potential of the form \x — y + QL\, where \x — y\ is the 
physical separation of the color charges, Q is an integer 
and L is the size of the box. These additional windings 
are due to the inability to gaug e away the gluon degrees 
of freedom in the problem 13, 0|. In effect, the treat- 
ment of Ref. [HI corresponds to including only the Q = 
sector. 

The study of these windings in Ref. [HI was motivated 
largely for theoretical reasons associated with QCD vol- 
ume reduction. Clearly in small boxes the two treatments 
might be expected to be quite different. A priori, how- 
ever, one might expect that the difference between cal- 
culations based on Refs. [r| and [l4| for a baryon mass 
would be quite small in cases where the box sizes are 
large enough to comfortably contain a baryon since the 
difference between the two is a finite volume effect. How- 
ever, as will be discussed in section llTII the differences can 
be quite large, even in such circumstances. Indeed, the 
differences are large enough that one might worry the dif- 
ferences suggest a possible problem with the formulation 
of Ref. [HI . Unfortunately, It is difficult to verify numer- 
ically that this is the case by going to exceptionally large 
box sizes. 

Given this situation, it is quite useful to have a regime 
in which one can independently compute the mass of a 
baryon to verify that the winding effects are as large as 
they appear to be for practical-sized boxes. The regime of 
heavy quark masses and large N c serves this purpose ex- 
tremely well. In this paper, we perform a highly accurate 
variational calculation of single baryons using Witten's 
mean-field argument [l7j . which is similar to an analo- 
gous 3+1 dimensional calculation that was done recently 
in Ref. [H)]. It was argued in Ref. 17 1 that baryons are 



composed of N c quarks in a color singlet state, where 
each quark moves in an average potential due to the re- 
maining N c — 1 quarks. In the non-relativistic limit with 
quark masses taken to infinity, the average potential has 
a simple form because the only contributions to the inter- 
action come from single gluon exchange. Note that the 
principal difference between the heavy quark system us- 
ing Witten's approach and the general result of Ref. 13 1 
is that the vacuum is trivial in the former case. 

We find that the single baryon results using Witten's 
variational formulation agree with the results based on 
Ref. [HI that incle the nontrivial effects of the winding 
with remarkable accuracy. These results support the ap- 
proach taken in Ref. [l3| that include a large effect of 
the windings, albeit only in the restricted regime of large 
quark masses. Our result is also quite useful in provid- 
ing a simple and physically intuitive picture for at least 
one regime of QCD in 1+1 dimension and, as such, nicely 
complements the Bringoltz method, which obtains results 
largely through numerical means. 

The study of the heavy quark regime in 1+1 dimension 
is useful for another purpose: the study of baryonic mat- 
ter at low density. Studies of low-density baryonic matter 
in the heavy quark regime of l arg e N c QCD were previ- 
ously done in 3+1 dimensions [18j |. The principal attrac- 
tion of this regime is that it appears to be tractable due 
to Witten's mean-field arguments and the non-relativistic 
nature of heavy quark dynamics. However, as discussed 
in Ref. [l8[ for the case of baryonic matter (as opposed to 
a single baryon) a plausible but ad hoc assumption was 
needed: that at low densities the ground state is com- 
posed of individual baryons each in a color-singlet state 
with the wave function of quarks for each baryon orthog- 
onal to the wave function of quarks in any other baryon. 
This structure appears quite plausible, but it is ad hoc in 
that one can construct mean-field configurations in which 
the entire state is color-singlet but cannot be broken into 
individual color singlet baryons. In 3+1 dimensions, it 
was not easy to verify whether this assumption was justi- 
fied. However, in 1+1 dimension, one can use the meth- 
ods of Ref. [l3[ to numerically compute the interaction 
energy and compare it with the result in the mean-field 
approach using the assumption of no hidden-color states. 
We find good agreement between them, which suggests 
that this assumption is indeed valid. This agreement in 
1+1 dimension greatly increases our confidence that the 
assumption made for the 3+1 dimensional case is valid. 

The paper is organized as follows. In the following 
section, we calculate the single baryon energy and wave 
function in the limit of large quark masses using Witten's 
mean-field argument. We compare this variational solu- 
tion to the numerical results based on the recent Hamilto- 
nian approach proposed of Ref. [l3[ and with a previous 
Hamiltonian in Ref. 141 that did not include the wind- 



ings. In the next section, we calculate the interaction en- 
ergy per unit cell of a low density baryon crystal, again 
using Witten's mean-field argument. We compare the 
interaction energies with numerical results for the many 



3 



baryon sector based on Rcf. [13J. The paper concludes 
with a brief discussion of the results. 



II. BARYONS IN THE COMBINED LARGE N c 
AND LARGE m q LIMITS 

A. The Energy Functional in the Mean Field 
Approximation 



The nonrclativistic mean-held energy functional for 
baryons in 3+1 dimensional large N c QCD for heavy 
quarks has been calculated in Ref. [l8[ based on Wit- 
ten's classic analysis in Ref. [17[. QCD in 1+1 dimension 
at large N c is generally simpler than in 3+1 dimensions. 
In 1+1 dimension, quarks are spinless, nonrclativistic 
fermions and gluon-gluon couplings can be eliminated by 
appropriate choice of gauge (e.g. the axial gauge as in 



refs. [8( and [1J]). Thus the leading order contribution 



to the energy comes from a coherent one-gluon exchange 
between pairs of -/V c quarks. Interactions involving mul- 
tiple gluons are suppressed at least at a relative order of 
Furthermore, quarks loops are suppressed at rela- 
tive order of at least jj- as a result of the large N c limit . 
However, for the case of heavy quarks, these simplifica- 
tions do not help. For heavy quarks, the interactions are 
pair wise color- Coulomb in any event so the lack of gluon- 
gluon interactions adds no new simplifications. Similarly, 
spin effects are suppressed for heavy quarks. Thus, the 
leading-order energy functional in the dual expansion, 
is derived in an essentially identical manner to the 3+1 
dimensional case [TJJ. However, the color-Coulomb po- 
tential in the context of one spatial dimension is linear in 
distance as opposed to dropping off as i . The strength of 

a 2 

the couplin g is chosen to be ^- , adopting the conventions 
in refs. 
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The energy of the baryon per color at leading order in 
4- and # is 



B. The Variational Calculation 

We choose to calculate the ground state baryon wave 
function i/> m i n (x) using a variational approach. From here 
on, ip(x) will be taken to be the ground state wave func- 
tion -0min(a;)- To proceed, we define a scaling variable R 
and a scaled wave function f(y) where y is the dimen- 
sionless variable defined as y = Then, the baryon 
wave function is: 



(3) 



The normalization of the wave function ip(x) is indepen- 
dent of the scaling variable R as long as f(y) is normal- 
ized, i.e. Jdy\f(y)\ 2 = 1. 

Next we define functionals for the kinetic and potential 
energies, T and V: 



T[f(y)] 



dy 



\d y f(y)\ r ' 



. . , \f(y)\ 2 \f(y')\\ ,, 
v[f{y)\ = / d y d y « w - y I 



(4) 



Rewriting the energy functional in terms of the new 
variables R, T[f] and V[/], we get: 



imw = M ba 



ryoi i 



(5) 



For simplicity, we choose a normalized functional form 
for f(y) to be a polynomial times a Gaussian: 



f(y) = 



1 + J2 a,iy 2i ) e 2 



N, 



<Jdy[ 1 + E a t y 2i 1 

i=l 



(0) 



(V>|#|V>) _ Af, 



Nc. 



d.r 



baryon 

N c ' 
\d x ^(x)\ 2 
2m a 



(1) 

, |#r)|W)| 2 , , 
dxdx \x — x I 



where m q is the quark mass, ip{x) is the normalized 
baryon wave function and A = g 2 N c . Minimizing the en- 
ergy functional with respect to the wave function while 
imposing the normalization leads to a Schrodinger equa- 
tion for the ground state baryon. 



= e 



t)' 2 

_ X 

2m q ' 4 



^ J dx'\i> min (x')\ 2 \x - x'\j V'min(x) 



(2) 



where e is the Lagrange multiplier that ensures the nor- 
malization of the wave function. 



where n indicates the number of terms kept. Only even 
powers of y appear in the wave function because the 
ground state is expected to be symmetric. We found 
that truncation at n = 9 was sufficient to obtain the en- 
ergy with a very high accuracy. We now have an energy 
functional that is a function of R and . The next step 
is to minimize the functional with respect to each of the 
variables, i.e. dh(ip\H\ip) = for he{R, a,i}. Minimizing 
with respect to R gives 



R 



mm} v 

Aro g V[/G0] 



(7) 



which yields: 



- N c (m q + \r^[f{y))vi[f{y)]m q h 2 ' 3 



(8) 
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To proceed we minimize (tp\H\ip) with respect to each 
ttj for i = 1 to 9 with ag chosen to be 1. Including R, 
there are a total of 10 parameters in the variational cal- 
culation. 

The binding energy of the baryon is the second term 
in Eq. (| 10[) and is characterized by the quantity T 3 Vs. 
It converges rapidly with increasing n, the number of 
parameters, flj, in the trial wave function. Minimizing 
this quantity yields 

oi w -0.0450660, a 2 « 0.0202162, 
a 3 ps -0.0021246, a 4 ps 0.0001965, 
a 5 ps -6.6218907 x 10~ 6 , a 6 ps -4.4275328 x 10~ 8 , 
a 7 ps 7.8179284 x 10" 9 , a 8 ps 4.2375112 x 10" 11 , 
-4.6069215 x 10~ 12 , 



(9) 



which fixes the wave function. 

It is sensible to present the results in units of A i.e. 
g 2 N c , which is the only scale present in large N c QCD. 
Hence masses will be written in the units of vA and 
lengths in the units of -7=. 

The ground state energy per color in units of vA is: 



0.25574 [ ^ 
,v A 



(10) 



The variational calculation converged well and the inter- 
action energy is estimated to be accurate to five signifi- 
cant figures. The parameter R is given by: 



n ( mm] y 

\Xm q V[f(y)]J 



the dimensionless Lagrange multiplier is 



1.0741901 



(11) 



(12) 



of the wave function - it is important to determine the 
correct behavior at asymptotic distances. 

To proceed, we study Eq. ([2]) in the limit where \x\ 3> 
\x'\. In doing so, we obtain a Schrddinger equation that 
captures the asymptotic behavior of the baryon wave 
function: 



3_ 

2m„ 



+ j\x\)^ asy (x) = eip aBy (x). 



The general solution is a linear combination of the 

Airy functions. Taking only the normalizable solution, 

i 

in units where x = (|-m 9 ) 3 x, wc find ip asy {x) = 
(|m 9 ) 6 ip asy {x) and e = (|) 



if e, we find that 



ipasyix) = aAi(23(|x| - e)) 



(14) 



where a is a constant. 

The asymptotic form is thus fixed entirely by the de- 
termination of the constant a. This can be done numer- 
ically provided that there is a region for which x is large 
enough so that the asymptotic form of the wave function 
of Eq. (fl"4"]l is accurate while simultaneously being small 
enough that the variational wave function is accurate. If 
so, one can match to the two forms in the region of si- 
multaneous validity of the two forms. It is easy to verify 
that such a region exists: a plot of the ratio of the two 
forms of the wave function must have a plateau with a 
value near unity over an extended region. Note that the 
size of such a region of stability depends on the accuracy 
of the variational wave function. It was largely for this 
reason, that we used as many terms as we did in the vari- 
ational wave function. Such a stable region was found in 
the vicinity of x — 3 and the asymptotic wave function 
takes the form: 



|[23(|S|- 



3 

r>]3 



^ ] - k 2^{\x\-e)]V* 
with k ps 1.21 and cps 1.07, 



(15) 



C. Asymptotic Behavior of the wave function 



where the asymptotic form of the Airy function has been 
used. 



The nature of variational calculations is such that they 
determine energies more accurately than the wave func- 
tions. The Gaussian form of the trial wave function we 
chose was very convenient and provided for an extremely 
accurate calculation for the energy. Moreover, with a 
sufficient number of terms, it provides a rather accurate 
form for the wave functions over the bulk of the region 
that contributes substantially to the energy. However, it 
does not capture the correct behavior of the wave func- 
tion at asymptotically large distances. At large x, the 
drop off of the true solution of Eq. ((2|) is slower than a 
Gaussian fl8j . Since one principal goal of this paper is to 
compute the interaction energy of low density baryonic 
matter - which is highly sensitive to the long range part 



III. COMPARISON OF SINGLE BARYON 
WITH MORE GENERAL APPROACHES 

In this section, we compare our single baryon energy 
results with numerical results based on the Hamiltonian 
approach of Refs. [l3| and [lij]. These approaches are 
more general than ours as they are valid for arbitrary 
quark masses. Our approach complements these; while 
the regime is far more limited, the approach is more 
transparent and far less reliant on numerical work. More- 
over, there is real value in having the ability to indepen- 
dently check the analysis of Ref. [l3| in a regime where 
the results of infinite box size are clearly under control. 
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Ref. [l3[ analytically finds the most general Hamilto- 
nian for large N c finite density QCD in the baryon sector 



J 



for arbitrary quark masses in a finite box. The Hamilto- 
nian for a box of size L is given by fl3j : 



8 M 



3 = 1 
M 



{ti:azd x p> (x,y) + c.c. + m q (x, y)} S(x - y) 



tr{p J (x,y)p j (y,x)} 



L + (x — y)(e m — 1J 



(16) 



2sin' 2 ( 



M 



r 



where the the trace is only over Dirac indices. The Hamil- 
tonian is given in terms of density matrices p> (a;, y) but 
can also be represented in terms of the conjugate space 
density matrix P®(x,y) as 



P j {x, y ) = J2 pQ ( x >y) e 

QeZ 



(17) 



The index Q represents integers and P®(x,y) is defined 
as: 

P Q (x 7 y) = ( (C[A]\^ (y)UHy, x + QL)^{x)\C[A\) ) A . 

(18) 

In Eq. (p~8|) . (O) a represents the average over all pos- 
sible gauge field configurations of some gauge-invariant 
quantity O 



(0) A = / d[A)0 , 



(19) 



U is a path-ordered (denoted V) spatial gauge field op- 
erator 



U(x,y)=Ve i9 ^ A{x ' )dx ' , 



(20) 



and |C[A]) represents appropriately normalized coherent 
states 

\C[A}) = e I dxdyC(x,y)g(x,y)^ _ ^ 

Note that |C[^4]) has implicit dependence on gauge fields 
(denoted A) via g(x,y), the gauge- invariant generators 
of the fermionic coherent group 1 1 31 ] 



g(x,y) = ^{x)U(x,y)ijj(y) 



(22) 



C(x, y) are the corresponding weights to g(x, y) and |0) is 
a color singlet reference state that satisfies the condition 
fln|0)=0j where a n is a member of a set of single particle 
annihilation operators chosen to ensure that the reference 
state has the baryon number of interest in the box. 

To recap, the diagonal projection onto coherent states 
defined in Eq. (|18[) is for a particular gauge configuration; 
however, reflects averaging over all possible gauge 
configurations. In this averaging, the only non-trivial 



structure that survives are the windings Q around the 
periodic box. 

Note that p^(y,x) has the following properties: 



p>{x,y) = p 3 {y,xY 
P 1 {x, y)p> (y, z)dy = p 3 (x, z) 



dx tr{fj (x,x) - p> vac {x,x)) = B. 



(23) 



Here p° vac {x,x) is the density matrix associated with the 
baryon Dirac sea and gives a formally divergent contri- 
bution to the integral in the absence of regularization. B 
is the baryon number in the box, which can be written 
in terms of b, the baryon density as B = bL. 

The baryon Hamiltonian of Eq. (|16p can only be solved 
numerically. To make the problem numerically tractable, 
we truncated the theory such that M in Eq. (|T6|) is kept 
finite. This truncation introduces corrections to observ- 
ables (including the baryon mass) that scale as jj . To 
obtain the observables in the full theory, the observables 
are calculated for different values of M and the result in 
the limit M — J- oo is obtained by extrapolation. 

To formulate the problem numerically, the Hamilto- 
nian is also discretized in space and the continuum results 
are obtained by extrapolation to the continuum limit. A 
latticized Hamiltonian using staggered fermions is dis- 
cussed in detail in Ref. |l3| . Note that the act of latti- 
cizing the theory introduces an ultraviolet cutoff thereby 
regularizing the theory and rendering the vacuum den- 
sity matrix, pl ac (x,x), finite. Here, we only present the 
form of the density matrix p J (x,y) latticized using stag- 
gered fermions. These matrices satisfy the properties in 
Eq. (|23p can be written in terms of a finite dimensional 
basis of orthogonal wave functions <^ (x) for each j with 
the indices n, x — 1,2...L S , where L s is the number of 
lattice sites. The discretized density matrix is 



P J (x,y) = <jy> n {x)(jy> n {yy 



(24) 



71=1 



with the following orthonormality constraints on the 
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wave functions: 



x=l 



(25) 



The 



sum from n = 1 to h^- in Eq. (f2"4")l is the lat- 
tice version of the Dirac trace of vacuum density matrix 
(f vac [x,x) in Eq. ([23]). The integral becomes regularized 
on the lattice and gives a finite contribution. 

The Hamiltonian of Eq. (fl6| was derived in the axial 
gauge representation of QCD It is well-known that 
in the Hamiltonian approach to either QCD or QED, 
Gauss' law appears as a constraint and not as an equa- 
tion of motion. Ref. [l5| explicitly constructed a QCD 
Hamiltonian, with the Gauss' law constraint built in a 
manner that only the physical states of the Hilbert space 
play a role. 

Ref. [l3[ , uses the constrained Hamiltonian of Ref. 
and the coherent states approach of Ref. [l6| to solve the 
Hamiltonian in the baryonic sector, which is the leading 
order physics in a expansion. This follows from the 
fact that the coherent states form an ovcrcomplctc set so 
that any physical state can be expressed as an integral 
over them and that the overlap between coherent states 
becomes sharply peaked as N c — > 00. 

The Hamiltonian has a number of interesting features. 
First, for translationally invariant states, the classical 
Hamiltonian for the baryons exhibits volume indepen- 
dence and partial volume independence [H, [lj| . Second, 
color charges interact via a linear Coulomb potential of 
the form \x — y + QL\, where L is the size of the box and 
Q refers to the number of windings around the box. The 
sign of Q determines the direction of the windings. The 
presence of these windings is a result of the inability to 
gauge away all the spatial gauge degrees of freedom in 
a finite box. Specifically, there are N c — 1 spatial gluon 
windings that remain in the axial gauge after the impo- 
sition of the Gauss' law constraint. 

Ref. |14| |. on the other hand, formulates the baryon 



Hamiltonian in infinite volume, where the windings rep- 
resented by Q are absent; the form of the Hamiltonian is 
exactly the same as that of Eq. (flBf except that Q = 0. 
However, in constructing the lattice regularized version 
of the resulting Hartree-Fock equations, Ref. 14| assumes 



that the correct approach is to use a box large enough 
such that finite volume effects on the baryon structure 
is small and with lattice spacing small enough that the 
baryon density can be easily resolved. For this purpose, 
Ref. [HI constructs a lattice potential, which has the form 
I a; — y\ in the continuum and infinite volume limits, while 
continuing to ignore the windings. Thus, the resulting 
latticized theory of Ref. [l4[ is equivalent to the M = 1 
truncation of the Hamiltonian in Eq. (|16[) and, for small 
boxes, may have significant errors. One expects these 
errors to go to zero as the box size goes to infinity. 



Comparison of Single Baryon Mass and Density 
Profile 



We find that in the limit of large quark masses, there 
is excellent agreement between the variational calcula- 
tion of section [TT] and a numerical calculation based on 
the latticized baryon Hamiltonian of Ref. In Figs. 

[1] and [2] we first plot the baryon mass minus the quark 
mass contribution (i.e. the interaction energy) as an ex- 
trapolation (in jj) for boxes of size WX = 2y/2n for 
= 100 and Ly/X = 1.6^/2^ for = 200. Note that 
the points plotted are the baryon masses in the contin- 
uum limit but for different finite values of M. We have 
points for M = 1, 5, 10, 15. 




FIG. 1: Baryon Mass Extrapolation for ^£ = 100 in a box of 
size Ly/X = 2v / 27r. Only interaction energies are shown. 



WcVT 



0.043 - 



0.042 - 



0.0 



1 
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FIG. 2: Baryon Mass Extrapolation for = 200 in a box of 
size L\f\ — 1.6\/2n. Only interaction energies are shown. 



The baryon masses (including numerical uncertainties) 
are ^= = 100.0550 ± 0.0001 for quark masses of ^ = 



100 and 



E 
N c \/X 



200.0440 ± 0.0003 for quark masses of 



-y= = 200. The uncertainties in the baryon masses are 
associated with taking the continuum and M — > 00 lim- 
its and also the systematic errors due to the numerical 
nature of the Hartrcc-Fock calculation. We find that the 
extrapolated baryon masses are in agreement with the 
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r~ i — m 

□ : M=15; /. s =170; La/ A =2a/2;t ; =100; 

Line : Baryon Density from Variational Method 




-2-10 1 2 



xa/T 

FIG. 3: Baryon Density for ^ = 100 with M = 15 and 
L a — 170 in a box of size Ly/\ = 2\/27r. 



r~ I — m 
□ : M=15; L s =170; La/ A =1.6yj2n ; =200; 

VT 

Line : Baryon Density from Variational Method 




-2-10 1 2 



xa/a~ 

FIG. 4: Baryon Density for ^ = 200 with M = 15 and 
L s = 170 in a box of size L^/ ^ X = 1.6a/27t. 



energies accurately for system with relatively large-sized 
boxes. From the plots in Figs. [3] and 0] it is clear that 
the baryons fit well within the boxes in the sense that 
the baryon density has dropped extraordinarily close to 
zero well before the edge of the box. Thus, one might be 
tempted to conclude that these boxes are big enough so 
that the effects of winding should be quite small. How- 
ever, this is not the case. Comparing the interaction 
energy of the M = 1 baryon mass (which we label by 
Eh) with the interaction energy after M — > oo extrapo- 
lation (labelled E H ) using the quantity i^~^ t ) , which 

measures the relative difference in the interaction pieces 
of the extrapolated and the M = 1 baryon masses. We 
also quote the error in the M = 1 baryon mass, which 
comes from systematic errors associated with the numer- 
ical nature of the Hartree-Fock calculation and the con- 
tinuum extrapolation. For a quark mass of ^= = 100, 

the M = 1 baryon mass is ^= = 100.05124 ± 0.00002 
and for quark masses of ^= = 200, the M = 1 baryon 
mass is —^7= = 200.04081 ± 0.00002. Thus, we find that 
the relative difference in interaction energy between the 
M = 1 result is approximately 7 per cent, which is quite 
significant given the degree to which the box fully con- 
tains the baryon. 

The large difference between the M = 1 and M = 00 
results show the practical importance of the windings in 
realistic calculations even for very large boxes. Of course, 
the effect of the winding is expected to be large for small 
boxes and was introduced largely to explain the behavior 
of the system in small boxes including large N c volume 
independence [Hj]. To some extent, the large effect of the 
windings on the boxes computed here is rather surprising. 
Indeed, they are large enough that one might question 
whether there is something wrong with the formalism in 
[l3j. However, the fact that the results extrapolated to 
M going to infinity agree so well with the results derived 
with the formalism of this paper, which are known to be 
valid for heavy quark masses, shows that this is not the 
case. 



masses from the variational calculation presented in Eq. 
(|10p . up to the numerical uncertainties of the calculation; 
the baryon masses from the variational calculation are 
= 100.055098 and = 200.043731 for quark 

masses of ^= = 100 and ^7= = 200 respectively. 

Next, in Figs. [3] and [U we plot baryon density pro- 
files for the two calculations with the line representing 
the result of the variational calculation and the dots rep- 
resenting the result of the numerical calculation of the 
general finite baryon density Hamiltonian. The densities 
are plotted for M — 15 and with the number of lattice 
points, L s = 170. The two results agree quite well in 
spite of the errors associated with discretization of space 
and truncation in the number of windings. 

It is useful to see the extent to which including the 
windings (i.e. M 7^ 1) is crucial for getting interaction 



IV. BARYONIC MATTER 
A. Introduction to the Multi-baryon Problem 

Now that we have established that Witten's heavy- 
quark formulation for single baryons is consistent with 
results based on the Hamiltonian approach of Ref . [l3| , 
we turn our attention to baryonic matter. It is believed 
that baryonic matter forms a crystal at large N c : the 
baryons are heavy and have a tendency lie in equilibrium 
at the potential wells created by the other baryons. 

We consider infinite nuclear matter at low enough den- 
sities such that only interbaryon interactions between 
nearest neighbors are relevant. 

Here, we assume that the many-baryon sector consists 
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of distinct color-singlct clusters of N c quarks with iden- 
tical wave functions, as in Witten's discussion of baryon 
scattering fl7l |. An analogous assumption was made in 
Ref. [lH for heavy quark, large N c QCD in 3+1 dimen- 
sions. This assumption is highly plausible but as noted in 
Ref. [HI, it has not been rigorously established. A possi- 
ble concern is that hidden color states could contribute at 
leading order. Thus, it is useful to test this assumption 
in a 1+1 dimensional setting. We can do that by compar- 



ing the results based on this assumption using our heavy 
quark formalism with numerical results using the meth- 
ods of Ref. [l3| , which are known to be valid at large N c 
generally and do not involve any additional assumptions. 

Let us begin by deriving the interaction energy per 
baryon in a 1+1 dimensional crystal at low density using 
the heavy quark non-relativistic formalism. The energy 
per unit-cell of the low-density baryon crystal in terms 
of a reference baryon is: 



E 

An 



dxdy|^(i)| 2 |^(y)| 2 |x-y| 



+ - / dxdy^i{x)^ i - 1 {y)^ i - l {x)^i{y)\x-y\ + - / dxdy^i(x)^ i+ i(y)^ i+1 (x)^i(y)\x - y\ , 



(26) 



where "J^ (x) is the wave function of the reference baryon, 
^i±i(%) are the wave functions of the nearest- neighbor 
baryons, the first term is the contribution to the intra- 
baryon kinetic energy, the second term is the contribu- 
tion to the intrabaryon potential energy, and the last two 
terms are the contributions to the interbaryon energy. In 
this expression both A and m q have been scaled out of the 
energy density using the following change of variables: 



A 



E= - 







rriq E . 



(27) 



By assumption, the density is sufficiently small so that 
the effect of interactions with next-to-nearest neighbors 
is negligible. 

The wave function for the baryon crystal is constrained 
by the Pauli exclusion principle, which precludes two 
fcrmions from simultaneously occupying the same po- 
sition. This implies that the wave function of a refer- 
ence baryon state ^>(x) must be completely orthogonal 
to the the wave function of the neighboring baryon states 
§(x + d) and $(x - d). 



= (1 + K)ipo(x) ~ ct(j)(x - d) - a<j>{x + d) 
with J dx\$(x)\ 2 = 1 

and / dx\^f {x)\\V(x±d)\ =0, 



(28) 



where 4>{x) is some normalized wave function. Here, we 
have assumed that densities are low enough that the ref- 
erence baryon wave function is dominated by the single 



baryon wave function ip(x) except for corrections cen- 
tered around x = ±d. The choice of k and a is deter- 
mined by the orthogonality constraint of the neighboring 
baryons and by normalization of the baryon wave func- 
tion. Wc find k and a in terms of 7 = J dxipo(x)4>(x), 
which characterizes the overlap at the same spatial points 
of the neighboring baryon wave functions and the quan- 
tity A' = J dx$o(x)<j>{x + d), which characterizes the 
orthogonality of the neighboring wave functions: 



a = — and k 

2 7 



AM 



where A 



A?_ 

4 7 2 

dx ipo(x)ipo{x ± d) . 



(29) 



We have disregarded 0(A 3 ) effects because A becomes 
exponentially small with decreasing density and in the 
low density regime, the higher order effects are paramet- 
rically small. Effectively, we are writing the interaction 
energy in the form A 2 f (where w is the characteristic 
width of a baryon and d is the separation of neighboring 
baryons) by developing an appropriate power expansion 
for /(f). 

To proceed, we first assume that <fi(x) = ipo(x), which 
appears to be the most efficient way to orthogonalizc 
neighboring crystal wave functions. We will show a pos- 
teriori that this assumption is correct in the sense that it 
will minimize the total energy up to corrections which are 
parametrically small at low density. The orthogonality 
constraint in doing so gives 



dx + d) = -A(d) + 0(A 3 



(30) 



As can be seen from the explicit form of A(d), which wc 
calculate in the following subsection, for parametrically 

-3 

low densities A{d) ~ e _d5 . Hence, it is plausible that 
in the low density regime, where only nearest-neighbor 
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interactions arc relevant, the choice of 4>(x) = ifto(x) is 
indeed the correct one. In the section HVDl we introduce 
corrections such that <j>(x) = rj(^ (x) + SA(x)), where A 
is a normalized wavefunction and argue that contribu- 
tions to energy at 0(5) do not contribute at the order to 
which we work. 



B. Calculation of overlap function A 

Before proceeding with the calculation of energy, it is 
useful to calculate the overlap function A, which is de- 
fined in Eq. <\'29\i. Most of the contribution to A comes 
from the localized baryon density at x = i. In the low 
density regime, the contributions to A, localized at x = 
and x = d are exponentially suppressed. The contribu- 
tion to the center comes from the tails of neighboring 
baryons, which have the form of Airy functions. Using 
the asymptotic form of Airy functions, we get 



the integral. This gives 



A= dx 



U2 g-i(2d-2g-22) 2 - I (-2e+2x)2 

2§7r(2d- 2e - 2x)3(-2e + 2x)i 



(31) 



where k rj 1.21 is the strength of the asymptotic wave 
function of Eq. (flip)) . 

The function in the exponent has a sharp minima at 
x = f and hence the integrand contributing to A has a 

maxima at x = ^. We use the standard steepest descent 
method to calculate the integral, which involves series 
expanding the exponent around the turning point, eval- 
uating the denominator of A at x = i and performing 



A(d) 



2§0F(J- 2e)"4 



(32) 



with a leading order correction term that is proportional 
The correction is suppressed for para- 



4(d-2€) 3 / 2 



t0 (d-2e) 7 /4 

metrically large distances, which is the relevant regime 
for this paper. 



C. Energy at low density 

Here, we calculate the interaction energy per unit cell 
of a low density baryon crystal up to G(A 2 ) and lead- 
ing order in O (f ) using units of Eq. ((211). Note that 
w is the characteristic width of an isolated baryon and 
d is the intcrbaryon separation. We divide the energy 
into two parts: the first line of Eq. (|26|) . which label the 
intrabaryon contributions, and the second line, which la- 
bel the interbaryon contribution. The intrabaryon con- 
tributions contain a contribution from the kinetic energy 
and another from a potential energy. Note that the in- 
trabaryon energy is not the same as in the free space 
baryon since the wave function has been modified due to 
the interactions with the other baryons. 

1. Interbaryon Energy 

The contribution to the interbaryon energy under the 
assumption that only nearest neighbor interactions are 
relevant is presented in Eq. (|33j) . 



^interbaryon 1 /* 1 

= - / dxdy %>(x)$(x + d)^>(y)ty(y + d)\x - y\ + - I dx dy * '(x)* (x - J)f (y - d)\x - y\ , (33) 



N c 2 J " y ' v ' x * " 2 

where ± d) = (1 + K) 2 ip (x)ip (x ± d) - all + k)^ (^)^o(») - «(1 + K,)tp (x ± d)tp Q (x ± d) 

--_ - - _ - - _ - ( 34 ) 

+ a'(il>o(x ± d)tp (x ± 2d) + i> (x)ip (x ± d) + ip (x)ip (x T d)) . 



The two terms in the equation correspond to the inter- 
baryon energies in a unit cell associated with two nearest 
neighbors of a reference baryon state ^f(x) in the crys- 
tal. Note that all terms in the interbaryon energy in 
Eq. Q33p are of the form J dx dy p(x)p(y)\x — y\, which 
is of the same structure as the Coloumb energy associ- 
ated with a one-dimensional charge density p{x). Thus, 
for example, the first term in Eq. ((3"3")) has an effective 
charge density ^(x)^(x + d). Moreover, note that in one 
spatial dimension, unlike in three spatial dimensions, the 
Coloumb energy of a charge distribution grows with the 
size of a system. 



In Eq. ((M)) . we present the effective charge density 
that determines the interbaryon energy in terms of the 
single baryon wavefunction if)o(x), a and n. Equation 
([29)) completely determines the two parameters. Since 
4> = tpo for a low-density baryon crystal, we have 7 = 1 
and A' = A. Therefore, a = y and K = jA 2 . 

To proceed, we note that the effective charge density 
relevant for the interbaryon contribution is localized. In 
Fig. [5j we show a schematic effective charge density plot 
characterizing the nature of ^(x)^(x — d) in Eq. (jM)) 
upto 0(A 2 ). Note that the densities are sharply peaked 
at x = , ^ , d. These peaks in Fig. [5] correspond to 
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FIG. 5: Schematic plot of the effective charge density 
'^(x)^i(x — d) associated with the interbaryon energy includ- 
ing the characteristic width of the relevant baryon densities 
in units of Eq. I)27p . The magnitude of the effective charge 
in each region is A. The total effective charges are A in the 
peak at x = % and — ^ each in the peaks at x = , d 



the first three terms of Eq. (fMf . The smaller peaks at 
x = and x = d correspond to the second and the third 
terms of Eq. (fM)) respectively. Since the isolated baryon 
wave funtion, ipa(x), is normalized, the terms each have 
a total baryon charge of — y. Furthermore, the width 
(in units of Eq. ([27]) ) of the peak at x = 0, which is 
proportional to the charge density, |-0o(a;)| 2 - is of 0(1) 
and does not scale with d as can be seen from the results 
of Eqs. ((6]) and (|11|) . The same properties are true 
for the charge density centered at x = d. However, the 
effective charge density centered around x = | behaves 
differently. For large interbaryon separation d, the width 
(in units of Eq. (f27|) ) scales as (d — 2e)*, which can be 
easily calculated from the standard steepest descent cal- 
culation in the integrand of Eq. (f2U)) . Also note that the 
total effective charge in the center peak is A for para- 
metrically small baryon densities. Although the width of 
the central peak depends on d (i.e. the baryon density), 
it is much smaller than the interparticle spacing at low 
density. Thus in the low-density regime, these three re- 
gions of charge are well separated from each other with 
a spacing much larger then the width of each lump. 

To calculate the leading contribution interbaryon en- 
ergy per unit cell of a baryon crystal in the low-density 
regime, the effective charge densities can be approxi- 
mated as delta functions with the appropriate baryon 
charges. The leading correction is the fact that the spa- 
tial extent of the charge density centered around x — | 
grows with d. The size of this correction can be calcu- 
lated using the approximation to the asymptotic form of 
the Airy function that was used to calculate the overlap 
integral A in Eq. (|33[) . This integral can be calculated 
easily using the method of steepest descent. 



\M*)\ 2 , 



the effective 



FIG. 6: Schematic plot of |#(x) 
charge density minus the single baryon charge density for < 
x < d including the total baryon charge in each region. 



For parametrically low baryon densities, the total in- 
terbaryon energy becomes: 



^interbaryon A^d ^^interbaryon 



Nr. 



(35) 



A TTiintcrbaryaii 

where N is the correction to the interbaryon en- 
ergy associated with the fact that the center peak scales 
with d: 



y\ ^interbaryon 

N~ c 



(d-2i)-iA 2 {d) 



(36) 



The interbaryon energy above has the form - charge 
squared times a width - which is as expected. Also, note 
that the energy contributions due to the finite extent of 
the two regions of effective charge at the the ends is O(d ) 
and hence is parametrically smaller at large d. 



2. Intrabaryon Energy: Potential Contribution 

In this subsection and the next, we calculate intra- 
baryon energy for the low-density baryon crystal at 
O (^) . There are two contributions to this, the potential 
and the kinetic. 

The potential contributions to the intrabaryon energy 
of the baryon crystal comes from Eq. (|37|) . The ef- 
fective charge density ^(x)! 2 is presented in Eq. (|38p . 
Additionally, we present in Fig. [5] a schematic plot of 
^(x)! 2 — |-0o(^)| 2 j the shift in the effective baryon charge 
density from the free case, ^(x)! 2 , minus the isolated 
baryon density, I^q^)! 2 , to emphasize the various charge 
enhancements. This quantity is even around x — and 
hence we only plot the region < x < d. 
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p intrabaryon 

potential = 1 / ^ ( _ } |2| _ _ _j > ^ 



W c 2 

where |*(x)| 2 = (1 + k) 2 ^ (^) 2 + a 2 (^o(i + J) 2 + ^ (^ - ^) 2 ) - 2(1 + /c)a(^ (z)V'o(x + <?) + $o{x)$o(x ~ <*)) 

2 - - - - (38) 
+ a (0o(x + d)ipo(x — a)) . 



To proceed, we calculate the effective charge density, 
^(x)! 2 , upto 0(A 2 ) to obtain the intrabaryon energy at 
the same order. Note that in Eq. (|38p . the first term, 
l^o (x)| 2 , of the effective charge density is centered around 
x = 0, the second term proportional to \tpo(x T <^)| 2 is 
centered around a; = ±ef, and the third term, which is 
proportional to ijjo(x)ijjo(x =F d), is centered around a; = 
±d. 

Furthermore, at 0(A 2 ) the intrabaryon energy contri- 
butions for the baryon crystal come from a cross-term be- 
tween the effective charge from the unperturbed baryon, 
■00 (x), and the effective charge due to the perturbation 
induced by the near- neighbor baryons. To the extent that 
d is large, the original charge is well-localized at the ori- 
gin, it creates an effective electric field which is constant 
for 3> 1 (or a potential linear in d). With a linear 
field dependence, the potential energy is sensitive to the 
center of charge the effective charge due to the perturba- 
tion. Thus, we can approximate the charge densities in 
Eq. (f38|) as localized delta functions with the appropri- 
ate strengths, which are associated with the total baryon 
charges of each of the terms. Now for the case of the first 
term, which has a strength of 1 + \A 2 and is localized 
at the origin, the delta function approximation is prob- 
lematic, since it is not in the region where the field from 
the unperturbed effective charge is linear. However, the 
potential energy from this term is of 0(A 2 w) : while the 
leading behavior scales as A 2 d. Thus, this term can be 



neglected as being parametrically small. The contribu- 
tions from the other two terms yield, 

^intrabaryon a2 J 

potential d /on\ 

~1T C — ■ (39) 

As noted above, there are corrections to this of relative 
order ^ from the perturbation in the effective charge due 
to k. There are other corrections due to the fact that 
unperturbed effective charge has a spread and hence the 
effective field is not precisely constant. These corrections 
are exponentially small in d. 



3. Intrabaryon Energy: Kinetic Contribution 

The kinetic energy contribution to the intrabaryon en- 
ergy at leading order in ^ is presented in Eq. (|40)) . Note 
that the first three terms are independent of the inter- 
baryon spacing d, whereas the remaining terms are not; 
in the low-density regime, the contribution from these 
terms independent of d are parametrically small. Fur- 
thermore, we ignore the 0(1) term because it is the con- 
tribution to the single baryon kinetic energy and not the 
crystal interaction energy. We also ignore the last term 
because it is not nearest neighbor. 



T^intrabaryon 1 /" 1 f 1 

= - / dx (d x V(x)) 2 = -(1 + nf / dx (d x Mx)) 2 + -a 2 / 



-a(l + K) J dx (d x ip (x))(d x ip (x + d)) - a{l + k) / dx (d x ip (x))(d x ip {x - <£)) 
+ a 2 dx (d x $ (x + d))(d x ipo(x ~ d)) 



(40) 



The remaining terms are trivially calculated using inte- integrand is sharply peak around x = ±|. 
gration by parts, the asymptotic form of the Airy func- 
tion in Eq. (TTJJ and the steepest descent method to eval- 
uate the integrals J dx (d x ^o(x))(d x 'iJjo(x=Fd)) 1 where the 
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The kinetic energy contribution to the interaction en- 
ergy then becomes 



E 



intrabaryon 
kinetic 

N r . 



A 2 d 



(41) 



4- The interaction energy 

Finally, putting all the interaction energies together, 
the interaction energy per unit cell of a low-density 
baryon crystal in the heavy quark mass limit up to 0(A 2 ) 
and to leading order in d is: 



E 



interaction 



Nr. 



pinterbaryon pnintrabaryon 
& J ^potential 

N~ r . W c 



~(d-2e)*A 2 (d) 

7T 



TTiintrabaryon 
kinetic 



(42) 



Equation contains as a factor (d-2lf/ i . We are 
interested in the limit of large d. So one could, in prin- 
ciple, expand (d — 2e)3 as times a series in powers of 
2- Moreover, since we are working in a regime, where d 
is parametrically large, all of the subleading terms in the 
series can be dropped. Hence, it may seem appropriate 
to simply use d* . However, it is actually prudent to keep 
all the terms in powers of 4. The reason is that there 

is a factor of (d — 2e)~3 coming from Eq. (|31[) . This 
means that our approximate expression for the interac- 
tion energy in Eq. (|4"2"j) breaks down as d — > 2e, yielding 
a diverging energy. This divergence goes as (d — 2e)~^ 
but if we replace (d — 21) i by d* in the numerator, it 
diverges much more slowly - as (d — 2e)~*. This slower 
rate of divergence suggests that the expression will re- 
main accurate for larger crystal densities. 



D. Self consistency 

Up until this point, we have assumed that the wave 
function for a baryon in unit cell was given by Eq. (|28|) 
with (f) = ipQ. Such a choice seemed natural in that it 
orthogonalizes the wave functions associated with differ- 
ent unit cells in a very efficient manner. It turns out 
that this ansatz can be shown to be correct: the correc- 
tions to it yield changes in the interaction energy that are 
parametrically smaller than the result given in Eq. (|4"2")l . 

To show this, we introduce corrections to 4>{x) of the 
form 



4>{x) =r)(Mx)+6A(x)), 



(43) 



where S is taken to be a small parameter, A(x) is an ar- 
bitrary wave function and 7/ is a normalization constant 
ensuring that <fi(x) is normalized. Without loss of gen- 
erality, one can always take take A(x) to be normalized 



and orthogonal to both ip(x) and tj;(x ± d). Any com- 
ponent in A(S), which are not orthogonal to tp(x) and 
if){x ± d) simply serves to renormalize r\ and n without 
altering the wave function. 

Our ansatz that <p = V'o is sclf-consistcntly correct if it 
minimizes the energy up to corrections, which are para- 
metrically small. This is true provided that for any choice 

of A, 



1 d ^i n t° rac tion 



d5 



O(A 2 d ) . 



(44) 



(5=0 



since the leading term in the energy is parametrically 
0(A 2 di). 

To proceed, one simply inserts the form in Eq. (|43[) 
into Eqs. (|28[) and (|26|) and expands to first order in 
S. Our self consistency condition will be satisfied if the 
terms of O(S) are parametrically of O(A 2 d ). There are 
a total of 31 terms in the interaction energy at order S 

- eight coming from the intrabaryon potential contribu- 
tion, 21 from the intcrbaryon energy and two from the 
interbaryon kinetic contribution. Of these, four terms 

- two from the intrabaryon potential contribution and 
the two from the intrabaryon kinetic contribution - are 
each of 0(A 1 ) 1 the other 27 terms are all of 0(A 2 ). It is 
straightforward to show that while each of the four terms 
of 0(A 1 ) is nonzero, their sum is exactly zero. 

Finally, there are 0(A 2 ) terms from the interbaryon 
and intrabaryon energies. It can be shown on a term- 
by-term basis that these terms all contribute at O(A 2 d ) 
or less. In doing so, the orthogonality of A(x)) with 
both 4>o{x) and V>o(S ± d) play a critical role since this 
implies that the integrated effective charge density due 
to cross-terms between A(\barx) and tpo(x) yield zero. In 
demonstrating that each term contributes at O(A 2 d ) or 
less, it is often necessary to separately consider the cases 
where A(x) is localized on a scale of w (that is, d°) or 
in which it is spread over a parametrically larger fraction 
of the unit cell since the reason that the contribution to 
the energy is O(A 2 d ) or less depends on the scale over 
which A(x) is spread. However, while the reason that 
a given contribution is parametrically small may depend 
on the scale of A(i), the fact that it is small does not. 

The fact that all contributions at O(S) either cancel 
or are parametrically of O(A 2 d ) or less, shows that the 
our ansatz was self-cosistcnt to the order at which we 
are working. Therefore, the energy per unit cell up to 
parametrically small correction is 



E 

N: 



Eo 

Nr. 



,Eh 

Nr. 



(45) 



where j^- is the leading order contribution, which is 
stated in Eq. (|4*2"j) and is the energy at 0{5 2 ). 
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V. COMPARISON WITH DIRECT 
NUMERICAL SOLUTIONS 

In this section, we compare the result of the low- 
density baryon crystal energy to the numerical results 
based on the Hamiltonian approach of Ref. [l3[, which 
as noted before is valid for large N c at any quark mass. 
Below, we plot the interaction energies for three differ- 
ent quark masses: ^= = 10,40, 100. Wc plot interaction 

V A 

energy (per color) versus the baryon density b = 
where B is the baryon number and L is the physical size 
of the box. The points represent the results from nu- 
merical finite density calculations based on Ref. [l3| ; the 
lines represent the results based on our heavy-quark-low- 
density formalism. 

The calculation was done for M = 10 and L s = 130 at 
a fixed volume for each of the quark masses. The 7 -^S= = 10 

vA 



calculation was done in a box of size L\f\ = 4^/2tt and 
the ^t= = 40, 100 calculations were performed in a box of 

size 

Obviously, there are numerical corrections to the in- 
teraction energy associated with discretization in space 
and the fact that M is finite. For the large quark masses 
involved the corrections associated with choosing a fi- 
nite M = 10 and a non-infinite number of lattice sites 
L s = 130 are relatively small The discretization 

correction depends on the number of lattice points that 
define the wave function of a single baryon. Since, we 
probed different densities in a box of the same size (for 
each quark mass) , the error associated with discretization 
is comparable for different baryon densities. 



2° 0.04 



FIG. 7: Interaction Energy with — 10, 40, 100 are repre- 
sented by triangle (dashed), square (dash-dotted) and circle 
(solid) respectively. The points are numerical calculations for 
the finite baryon density 't Hooft model using the method of 
Ref. [U and the lines are the result from the heavy-quark-low- 
density treatment developed in this work, b = -^j^ (where B 
is the baryon number in a box, and L its length) is a dimen- 
sionless measure of baryon density. 

We find that the agreement between the heavy-quark, 
low- density calculation of Eq. (|4"2"j) and the results based 
on Ref. [13| is remarkabably good for small b. Indeed, 



it is surprising just how well it works. In deriving Eq. 
(|42p. rather strong assumptions were made - including 
the fact that contributions of O(A 2 d ) were neglected 
as compared to terms of 0(A 2 d*). One might expect 
these to become very important as b increases. Note from 
Fig. [3 that the approximation breaks down dramatically 
at certain fixed values, which depends on the quark mass. 
As noted above, this break down occurs at d = 2e. It is 
remarkable, however, how accurate the large-mass-low- 
density expression remains even quite close to the break 
down value. 



VI. DISCUSSION 

In this paper, the nonzero baryon number sector of 
large N c QCD in 1 + 1 dimension in the limit of large 
quark masses was explored. The energy and density pro- 
file of individual baryons for different quark masses and 
the interaction energies of low density baryon crystals 
were computed. These obscrvablcs were calculated di- 
rectly using Witten's mean-field approximation, which 
simplifies greatly for infinite quark masses [TtJ since the 
result depends entirely on interactions involving single 
gluon exchanges. This approach, although valid in a lim- 
ited regime of quark masses has the virtue of providing a 
simple intuitive picture for nonzero baryon density which 
complements the numerical approach Ref. fl3| . The com- 
puted mean-field energy and baryon density profiles com- 
pare very well with ones calculated numerically based on 
the Hamiltonian of Ref. [l3[ . The agreement in the large 
quark mass regime is excellent, and demonstrate the sig- 
nificance of the effect of all the gluon windings even in 
cases where the baryon density is exceptionally small at 
the boxes' edge. 

The analysis of low density baryonic matter developed 
in this paper relied on an ad hoc assumption that the 
hidden-color states - in which the wave function cannot 
be factorized into the product of color singlet baryons, 
play a subleading role at large N c . This assumption is 
highly plausible; however, unlike for the single baryon 
case, has not been explicitly demonstrated to be correct 
[ijj. The fact that it agrees so well with the numerical 
solution, which was derived to be valid at large iV c with- 
out making such an ad hoc assumption, is strong evidence 
that the assumption is in fact valid. Analogous reasoning 
was used for the case of large N c QCD with heavy quarks 
in 3+1 dimensions [l8| in the low density regime; the suc- 
cess of the approach in 1+1 lends credence to the 3+1 
calculation for which there is no numerical cross-check. 

Finally, the success of our low density approximation 
with its ad hoc assumption in matching the numerical 
results based on Ref. [13| , which was designed to work at 
large N c without such ad hoc assumptions, suggests a way 
forward in formally establishing the validity of the ap- 
proximation used here and quite conceivably for the 3+1 
dimensional case. Note that the approach of Ref. [l3| 
was based on properties of a class of coherent states. It 
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is highly plausible that a formal treatment based on the 
coherent states approach of Ref. [l3| but taken for the 
heavy quark limit will yield the formalism developed here 
and thus replace our ad hoc assumption with something 
more solid on theoretical grounds. Moreover, it is plausi- 
ble that some extension of such a method may work for 
the 3+1 dimensional case as well. While the general 3+1 
dimensional case has a very complicated structure even 
at large N c due to the gluodynamics, in the heavy quark 
limit the gluons only contribute to produce static color 
Coulomb interactions and one might imagine an appro- 



priate coherent state formalism applying for the quark 
fields. These possibilities will be explored in future work. 
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